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Simple bounds are established for the solutions of second-order homogeneous linear difference 
equations in ranges in which the solutions are exponential in character. The results are applied to a 
recent algorithm for the computation of subdominant solutions of second-order linear difference equa- 
tions, homogeneous or otherwise. Strict and extremely realistic bounds are obtained for the truncation 
error associated with the algorithm in a number of examples, including Anger-Weber functions, Struve 
functions, and the solution of a differential equation in Chebyshev series. 
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1. Introduction 

Consider the difference equation 

Pr+1 _ bp r + Pr-\=0, 



(1.01) 



in which b is a real constant exceeding 2. The solu- 
tion may be expressed as 

= Pi-^zPo xr , XiPo-pi xr 



Ai A? Ai A2 



K. 



where po and p\ are at our disposal, and 



^=' 2 b + 



VF^ 



^ = - 2 b- 



4 



1. 



Eventually \p r \ grows in proportion to Ai as r— * °°, un- 
less p\ = X 2 po- The larger the value of 6, the faster is 
the rate of growth. 

Now consider the equation 



Pr+l — 6rPr + Pr-l = 0, 



(1.02) 



in which b\, b 2 , 63, ... . is a given sequence of real 
numbers, not necessarily monotonic, each of which 
exceeds the b of eq (1.01). If the same values of po and 
Pi are used for the two difference equations, then it 
is reasonable to expect — and not difficult to prove — 
that the solution of (1.02) generally increases at a 



faster rate than that of (1.01). Similarly, if, for every 
r, b r < B, where B is another constant, then we ex- 
pect the solution of (1.02) to grow at a slower rate than 
the solution of 



p r+ i — Bp r 4- p r -i = 



(1.03) 



having the same initial values. 

The first purpose of the present paper is to generalize 
these intuitive results and present them as precise 
theorems. We shall allow the coefficient of p r -\ in 
(1.02) to take values other than unity, and we shall 
not insist that the solutions of difference equations 
under comparison agree at both r=0 and r=l. Some 
related results for the tails of continued fractions have 
been given by Blanch [1] l and Merkes [2]. 

The second purpose of this paper is to apply these 
theorems to a recent algorithm for computing sub- 
dominant solutions of homogeneous or inhomogeneous 
second-order difference equations [3], with a view to 
obtaining bounds for the truncation error. 



2. Comparison Theorem 



Theorem 1. Let 




qr+i = /3rqr — « r q r -i, 


(2.01) 


Qr+i = B r Q r — A r Q r _i , 


(2.02) 


here {a r }, {ft}, {A r }, {B r }, r=l, 2, . 


., are se- 



1 Figures in brackets indicate the literature references at the end of this paper. 



quences of numbers such that 
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A r ^a r ^O, B r -A r ^/3 r 



1. 



// 



Qi — qi ^ Qo — qo ^ 0, and qi ^ max (q , 0), 

then q r and Q r are nondecreasing functions o/r, and 
q r 2=Q r ,r=0, 1, .... 
To establish this result, we begin with the identity 



(2.03) 


respectively, and 






a— inf a r , 


b = a-hinf (b r — a r ), 


(2.04) 


A=sup a r , 


B = A 4- sup (b r — a. 



(3.06) 



(3.07) 



Before entering into the proof of this theorem, we 
observe that the numbers b and a always exist and 
are subject to the inequalities 



q r +\ — q r 



(fir — ctr—l) qr + ot r (q r — q r -\ ) 



(2.05) 



6^1 + a^l. 



(3.08) 



Moreover, 



and observe that both /3 r — ot r — 1 and a r are non- 
negative. From (2.04) we see that q x ^ and q x — qo ^ 0. 
Hence from (2.05), with r= 1, it follows that q 2 — qi^ 0. 
Continuing this argument by induction, we see that 
q r +\ ^qr, r=0, 1, . . . . Similarly from the identity 

Qr + i-Qr=(B r -A r -l)Q r ^A r (Q r -Q r - l \ (2.06) 

we deduce that {Q r } is a nondecreasing sequence. 
Next, subtraction of (2.05) and (2.06) yields 

(Qr+i — q r +i)—(Qr — qr)= (j3 r — a r — 1) (Q r — q r ) 

+ {(Br-Ar)-~( r 3r-Ctr)}Qr 

+ O r {(Q r -qr)-(Qr-l-qr-l)} 

+ (A r -a r )(Qr-Qr-l). 

That p r is a nondecreasing function of r is an im- 

Assume that r ^ 1 and Q r - q r ^Qr-\-qr-\ ^0-as mediate consequence of (3.02), (3.03), and Theorem 1. 

is certainly the case when r=l in consequence of We shall establish the inequalities (3.04) in turn by 

(2.04). Then using (2.03) and the facts that Q r ^ and further use of this theorem. 

Q r -Q r _ x ^0, we deduce that (Q r +x-q r +\)-(Qr-qr) Suppose first that k ^ X. Define q r = pok r . Then 

^0. Hence by induction Q r ^ q r . This completes the <7o = Po, 9i=PoX and, in consequence of (3.05), 
proof. 



X = 2&+ -J j b 2 — a^ -a + --\-\-a— ~ | = max (a, 1). 

(3.09) 

On the other hand, either or both of the sequences 
{a r } and {b r — a r } may be unbounded, causing either 
B or both B and A to be infinite: in this event A is unde- 
fined and the second of the inequalities (3.04) is inap- 
plicable. Except in these cases, we have 2 



A^a, B^b, A^A, 



A^X. 



(3.10) 



3. Bounds for the Solutions When the a r 
Are Nonnegative 



q r +i = bq r -aqr-i (r^ 1). 



(3.11) 



Theorem 2. Let 

p r +i =b r p r — a r p r _i. 

where {a r } and {b r }, r=l, 2, . . ., 
numbers such that 



(3.01) 



We apply Theorem 1, with (3.11) filling the role of 
(2.01), and (3.01) that of (2.02). Conditions (2.03) are 
satisfied in consequence of the definitions of a and 6, 
and (2.04) is satisfied since /c^X^l. Therefore 
p r = q r — PoX r , in agreement with the first of (3.04). 
Now consider the case when k < X. Define q r = poK r . 



sequences of Then *> = "^ * = ^' and 

q r +i = bq r — dqr-i (r^ 1), 



(3.12) 



b r ^l + a r ^l, (3.02) 

and suppose that 

p > 0, k = pJpo ^ 1. 

Then p r is a nondecreasing function ofr, and 

p {min (k, X)} r ^ p r = Po{max (k, A)} r 

(r = 0, 1, . . .). (3.04) 

Here X, A are the largest roots of the equations 

X 2 -bX+a = 0, A 2 -BA + A=0, (3.05) 



where a = min (a, k) and b = k-\-clk~ 1 . The desired re- 
sult p r = q r will follow if Theorem 1 can be applied 
to eqs (3.12) and (3.01). Conditions (2.04) are clearly 
(3 03) satisfied; the remaining conditions are 



5^0, 



ar^b-a^l. 



The first two inequalities follow from the definitions 
of a and a. And since 

6-a-l-(/c-l)(l-5K" 1 ) ^0, b-a^br-a r , 



2 The last of (3.10) may be verified by expressing 2\ = b+ (b 2 — 46 + 4c) 1/2 , where 
= inf (6r"~«r)S 1, showing that dk/db and dk/dc are both positive. 
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it will suffice to show that 

b — a ^ b — a. 



(3.13) 



Theorem 3. Let 

p r +i = b r p r +a r p r _,, 



(4.01) 



If d=K, then b = K-\-l and (3.13) follows immediately where 
from (3.08). If, on the other hand, a = a< k, then 



[b-a)-{b-~a) = b- 



clk 



K + - 



(X-„)(l-jL)>0. 



since a < k and 1 ^= k < X. This completes the proof 
of the left-hand inequality in (3.04). 

To establish the right-hand inequality in (3.04), sup- 
pose first that k^A. Defining Q r = PoA r , we have respectively, and 
Qo = Po, (>i=poA, and 



a r ^0, b r ^0, p >0, * = pi/po£l. (4.02) 

Then 

p {min (k, X)} r ^ p r ^ p {max (k, A)} r 

(r=0, 1,. . .), (4.03) 
where X, A are the largest roots of the equations 

X 2 -b\-a = 0, A 2 -BA-A = 0, (4.04) 



Q r+x = BQ r -AQ r - 



(rfcl). 



(3.14) 



It is readily verified that Theorem 1 is applicable to 
eqs (3.01) and (3.14). Hence p r = poA r , in agreement 
with (3.04). 

Now consider the case when k > A. Define Q r = PoK r , 
so that Qo = po, Q\=Pu an d 



a = inf a r 



A = sup a r , 



6 = inf b r , 
B= sup 6 r 



(4.05) 
(4.06) 



To prove the left-hand inequality in (4.03) when 
k i^ X, assume that p r ^ pnX r and p r -\ = poX r_1 , — as 
is the case when r=l. Then from (4.01) we have 

Pr+i ^ bp<s\ r + apoX r ~ 1= poX r+1 , 

as required, 
where #=k + ^k _1 . The desired result p r = Q r will Now consider the case when k < X. The roots of 

follow immediately if Theorem 1 is applicable to (3.01) the quadratic equation x 2 — bx — a = are X and — a/X. 
and (3.15). Conditions (2.04) are obviously satisfied. Since —a/X < k < X, it follows that 
Conditions (2.03) demand that 



■■BQr-AQr-i (rSl), 



(3.15) 



< 2 -OK-a<0. 



(4.07) 



^ ^a r ^0, B-A^br-ar^ 1. 

Assume that p r = poK r and p r _j ^ poK r ~ l — as is the 
In consequence of (3.02) and (3^07), these inequalities case when r=l. Then from (4.01), (4.05), and (4.07), 



will follow if we can show that B ^ B. Now 



deri\ 



B-B=k + - 



+£]=(fC-A) ( 



Ak 



Pr+l ^ bp K r + ap K r ~ l > p K r 



This is nonnegative because k > A, k ^ 1, and A i^; 
compare (3.10). The proof of Theorem 2 is now 
complete. 



4. Bounds for the Solutions When the a r 
Are Nonpositive 



If the quantities a r are negative then Theorem 2 is 
inapplicable. The analysis of this case is somewhat 
easier, however, and Theorem 1 is not needed. For 
convenience, we replace a r by — a r . 



The right-hand inequality in (4.03) may be estab- 
lished in a similar way. 



5. Anger-Weber Functions 

In the remaining part of this paper we show how the 
preceding theorems may be applied to the truncation 
errors associated with the algorithm of [3]. 

Consider first the computation of the Anger-Weber 

function E,(l) given in [3], section 6, Example 1. The 
truncation error of the approximation y\} 4) is given by 3 



3 [3], section 5. 
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4 14 >^E,.(l)-y r 14) =£i4/v, 



where 



£14=2 



e s 



(5.01) 8.248 X 10~ 17 by direct numerical summation of the 

terms in the expansion (5.02). Our strict bound (5.08) 

exceeds this estimate by about 7 percent, which is 

very satisfactory in view of the somewhat crude step 

/r ^2^ leading from (5.05) to (5.06). 

Strict bounds for the errors € ( r 14) (r^l4) may be 
computed by substituting the numerical value (5.08) 
p r satisfies (3.01) above with a r =l, b r =2r, p = 0, for E u in (5.01) and using the values of p r given in [3], 
Pi = l, and table 1. Obviously they, too, will overestimate the 

„ . ,_ ~ ox actual truncation errors by about 7 percent. 

e r = er-l - drPr, d r = - (2/ir){ 1 - (- l) r } , (5.03) 



with eQ= Eo (1). 



6. Struve Functions 



We apply Theorem 2 of section 3 to the sequence As a second example, consider the Struve function 

{pr+14}, using the numerical entries given in [3], H r (0.1) computed in [3], section 6, Example 2. The 

table 1. Working to five significant figures (for the sake relevant formulas in this case are 
of illustration), we have 



K = Pl5 / Pl4 = 1.1127 X 10 15 /(3.9793 X 10 13 ) = 27.962, 
a =l, & = l + inf (2r-l)=30. 

Hence 

\ 2 -30\ + l = 0, \ = 29.967>k. 
Therefore 

p r ^Pi 4 K r - 14 (ri= 14). (5.04) 

Next, from (5.03) we have 

e s =e 14 +(4/7r)(p 15 + Pi7+. . .+P.0 (s > 14) , 

(5.05) 

where s' = s if s is odd, and 5' = 5—1 if s is even. 
Since p s is a nondecreasing function of 5 (Theorem 2), 
we have 



€< r 15 >=H r (0.1)-yJI 5 > = i? 1 5Pr, 



€ S 



PsPs+1 



(6.01) 



and 



i = 20rp r — p r _i, 



-e r -i — 



(0.05)' 



*H) 



V^r 



(6.02) 



Applying Theorem 2 to the sequence {p r +is} using 
the numerical entries in [3], table 2, we have 



K = p 16 /Pi 5 = 299.99643. . ., 
o=l, 6=l + inf (20r-l): 

X=160-hV25599>K. 



320, 



Hence 



e s ^ e 14 + (4/tt) (s-U)p s (* B 14) . 
Application of this inequality to (5.02) gives 



(5.06) 



s ff 4 PsPs+i rr \pi6 Pi 7 Pis / 



p r ^p 15 /c r - 15 (rS15). 
Inspection of table 2 of [3] suggests that 
|e r |^e 15 (rS15). 



(6.03) 



(6.04) 



Substituting (5.04) on the right-hand side of (5.07), 
summing the resulting series and using numerical 
entries given in [3], table 1, we obtain 



Deferring the proof of this inequality for the moment 
(5.07) and substituting (6.03) and (6.04) in the second of 
(6.01), we derive 



l*i. 



Cl5 



P?5 (k-K' 1 ) 



— = (0.23277 966 . . .) X lO" 62 . 

(6.05) 



E 14 ^ 



e 14 



+ - 



P 2 4 (k—K' 1 ) 7Tpu(K—l 



- = 8.824 XlO" 17 . 



The actual value of £15 may be estimated by compu- 
tation of P17, pig, . . . and direct summation of the 
(5.08) expansion in (6.01). This gives (0.23277 943 . . .) 
X 10~ 62 . Therefore strict bounds for e ( r 15) (r^ 15) com- 
The actual value of Eu was estimated in [3] to be puted from (6.05) and the first of (6.01) will exceed 
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the estimated errors by less than one part in 10 6 , a in table 4 of [3] for r— 0, 1, . . . , 8. Applying Theorem 
striking example of the power of Theorem 2. 
It remains to establish (6.04). Write 



ttr =(0.05)W{V^r(r+|U 



2 to the sequence {/v + 7}, we have 

K = p s lp 7 = 5.4389, 
correct to five significant figures, and 



(7.04) 



that 



e r — ei 5 — ^16 — ^17 — 



u r (r^l6). 



The p r and u r are all positive, in consequence of (6.03). 
From the first of (6.02) we see that 



p r+ i < 20rp r . 

Accordingly, we derive 

a r+1 <2ra r /(2r+3), 

and thence 

/30 30 32 , 30 32 34 



. r 2r-l 15 . 15 , . r 10r+l 96 

r ^s2r-\-l 17 17 r ^8 2r+l 17 

Accordingly, 

17\ 2 - 96X + 15 = 0, \ = 5.4862 > k. 
Hence 



Pr^p 7 K r - 7 



(rB7). 



(7.05) 



Substitution of the last inequality in the second of 
(7.03) yields 



2 Ur K " 15 133 + 33 '35 + 33' 35 '37 + 



• - -)= 30 ^ E7 " ^(li + T^+T^+ • • •) 



< T^T7 777=0-67713 X10- 10 . 

15p^ (k — k *) (7.06) 

The corresponding bounds for the errors c|) ( r 7) are given 
in the second column of table A below, in units of the 
tenth decimal place. They exceed only slightly the 
estimates for the errors obtained by subtracting the 
entries in the columns headed f[ l2) and / ( r 7) in table 5 
of [3], and given in the third column of table A. 



Table A 



the last relation being a consequence of Gauss' 
formula for the hypergeometric function F (15, 1; 
16* ; 1); see [4], eq (15.1.20). Since 

30ai5 = 0.01418 23625 . . . 

the stated result now follows. 

7. Clenshaw's Numerical Method for Ordinary 
Differential Equations 

As a final illustration, we consider a recurrence 
relation arising in the solution of a differential equation 
by the method of Clenshaw; compare Example 4 of 
[3], section 11. Suppose first that we are interested in 
the solution of the equation 

(2r-l)/r-i-12r/ r +(2r+l)/ r+ i = (7.01) 

satisfying the conditions /o=l and/ r — »0 as r— »<». 

Approximations to f r are provided by the function 

fi 7) given to six decimal places in table 4 of [3] and to c „ , . , . . 

J r , . , , . 1 1 r r 1 r Secondly, suppose that we are interested in the 

ten decimal places in table 5 of the same reference. golution of (m) which satisfies {he condition 

Let us seek bounds lor the magnitude ol the truncation 



r 


0.67713p r 


iQ'o^m 
(estimated) 


0.11303/* 


-3735//' 


10V r 7 > 
(estimated) 














-3735 


-3655 


1 


1 


1 





-322 


-315 


2 


3 


3 





-41 


-40 


3 


13 


13 


2 


-6 


-3 


4 


63 


63 


10 


-1 


9 


5 


325 


324 


54 





54 


6 


1723 


1715 


288 





286 


7 


9268 


9228 


1547 





1540 



error 

<K. 7) =/r-/< 7 >. 

As in the other examples, we have 



2yo + yi + y2 + y3 + . 



= 1. 



(7.07) 



\j.\jZ) (This is the form in which the example was originally 
proposed in [3].) Then y r is related to/ r , defined above, 
by 



♦«-**, ft-Sj^r 



(7.03) 



In the present case p r satisfies the difference eq (7.01), 
and e r = l/(2r-h 1). Numerical values of p r are given 



y r =f r /F, F=*/ +/i+/,+ . . . . (7.08) 
We seek bounds for the errors 

6 (7) = Jr _y r 7) ( r =0, 1, ... 7), (7.09) 



j 280-2 10 0-68— 2 
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where 



cr 7 < 13403 X 10- 



yip=f(i>IF 7 , F 7 =l/< 7 >4^ 7) +/ ( 2 7) + - . .+/£ 



Using (7.02) and (7.08), we find that 



A7) = A7,l) — A7,2) 



where 



and 



,(7) 



: (7, l) = Vr 

F 7 + o- 7 ' 



A7,2) = 

*- r 



gt£g 

F 7 (F 7 + or 7 )' 



(7.10) 



(7.11) 



(7.12) 



(7 7 = H!) 7) +^ 7) -+-^ 7) +- • - + <H 7) + 

+/8+/9+/10+. • . . (7.13) 

In order to assess an upper bound for cr 7 , we need 
an upper bound for 



fr = Pr £ 



1 



r (25 + VjPsPs + : 



(7.14) 



when r^8; compare [3], eq (5.03). To construct this, 
we make a second application of Theorem 2 to the 
sequence {pr + 7}. We have 

a 2r-l 10r+l 

^ = SU P ^TTT = ! ' B = l + sup 9 r -ui = 6 - 

r ^ 8 2r-h 1 r^ 8 ^r-r 1 

Accordingly 

A 2 -6A + 1 = 0, A = 3 + V8 = 5.8284>k, 
where k is defined by (7.04). Hence 

Pr ^p 7 A r - 7 (rS7). (7.15) 

Substituting this result and (7.05) in (7.14), we find that 



A 



/-T-2 



1 



p 7 ^ r (25"+l)^ 



< 



and thence 



A' 



p 7 (2r J rl)(K-K- 1 )K 2r - 



(rS7), 



-/ r< 17p 7 (/<-/<- 



1 )(k 2 -A) 



= 2007X10- 



correct to ten decimal places. Adding this value to 
10" 10 times the sum of the entries in the column of 
table A headed 0.67713p r , we derive 



compare (7.13). 

Since F 7 and cr 7 are both positive, we obtain from 
(7.03), (7.06), (7.12), and the value F 7 = 0.59907 ex- 
tracted from [3], table 5 



P (7, 1) : 



F 7 



E 7 p r 

F 7 



^ 0.11303 Xl0~ 9 pr 



and 



47, 2) ^ £lf(i) ^ 3735 x l0-9/ ( r 7) , 



both e ( r 7, 1} and e ( r 7, 2) obviously being nonnegative. The 
values of 0.11303p r and — 3735/ ( r 7) are given in the 
fourth and fifth columns of table A. Comparison with 
the estimated values of 10 9 e ( r 7) , extracted from table 5 
of [3] and given in the final column, is again fully 
satisfactory. 

8. Summary 

In the first part of this paper (sees. 1-4) simple 
lower and upper bounds are established for the solu- 
tions of second-order homogeneous linear difference 
equations, in ranges in which the solutions are ex- 
ponential in character. 

In the second part (sees. 5-7), the results are applied 
to an algorithm for the computation of subdominant 
solutions of second-order difference equations which 
was introduced recently by the writer. It is shown by 
means of examples how to compute strict and ex- 
tremely realistic bounds for the truncation error 
associated with the algorithm. The bounds depend 
only on simple properties of the coefficients in the 
difference equation, and are independent of asymptotic 
theories of the solutions. 



The writer is grateful to C. W. Clenshaw and G. F. 
Miller for valuable criticisms of an earlier form of 
this paper. 
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